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Abstract 

Superconducting and antiferroquadrupolar (AFQ) transitions in a Pr-based compound 
PrRri2Zn2o have been found to occur simultaneously at T c =Tq=0.06 K. The superconducting 
transition manifests itself by zero resistance and large diamagnetic susceptibility. The specific heat 
exhibits a Schottky anomaly peaking at 14 K and magnetization curves measured at 2 K show 
anisotropic behaviors. The analysis of these data indicates that the crystalline electric field (CEF) 
ground state of the trivalent Pr ion is the non-Kramers T3 doublet with the quadrupolar degrees 
of freedom. A sharp peak in the specific heat at 0.06 K has been attributed not to the supercon- 
ducting transition but to the AFQ transition because the ordering temperature Tq decreases in 
B || [100] but increases in B\\ [110] and B || [111] with increasing B up to 6 T. This anisotropic 
behavior of Tq(B) can be well explained by a two-sublattice mean-field calculation, which corrob- 
orates the AFQ ordered state below Tq. The entropy release at Tq is only 10% of Bhx2 expected 
for the T3 doublet, suggesting possible interplay between the quadrupolar degrees of freedom and 
the superconductivity. 

PACS numbers: 71.70.Ch, 74.70.Dd, 75.20.Hr, 75.25.Dk, 75.30.Kz 
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I. INTRODUCTION 



Rare-earth intermetallic compounds have attracted much attention because of a variety 
of phenomena originating from 4/ electrons such as heavy-fermion state, superconductivity, 
Kondo effect, and multipole ordering. In these systems, the total angular momentum J of 
the 4/ electrons is the good quantum number because of the strong spin-orbital interaction. 
Therefore, the observable is not the orbital degrees of freedom but quadrupoles which are 
described as second-order tensors of J. When the dipole magnetic moments are quenched 
but the quadrupoles remain active under the crystalline electric field (CEF), which is pos- 
sible in non-Kramers ions at a cubic point group local symmetry, the quadrupoles often 
play an important role in forming exotic electronic ground states such as antiferroquadrupo- 
lar (AFQ) state with a staggered quadrupolar component Ql and a non- Fermi liquid (NFL) 
state attributed to two-channel (quadrupole) Kondo effect. 2|-|4j Furthermore, the feasibility 
of superconductivity mediated by quadrupolar fluctuations in the cubic Pr-based supercon- 



ductor PrOs4S 
measurements. 



3i2 with T c =1.5 K has been pointed out by neutron scattering and NQR 
5143] Another well-known cubic system is PrPb 3 which undergoes an AFQ 
transition at Tq=0.4 K. [g, Q| In this system, the CEF splits the ninefold multiplet of J=4 
into four multiplets. The CEF ground state of T3 doublet carries no magnetic dipole moment 
but electric quadrupole moments 2 = (3Jf — «7 2 )/2 and 0|=a/3(J^ — Jy)/2. The former 
has been found to be the AFQ order parameter by the combined analysis of magnetization 
and neutron diffraction experiments. IQj, 111] The quadrupole moments are aligned with an 
incommensurate sinusoidally modulated structure even in the ground state. Thereby, the in- 
direct Ruderman-Kittel-Kasuya-Yosida-type interaction between the quadrupoles plays the 
essential role. On the other hand, substituting La for Pr in PrPb3 destroys the AFQ or- 
der at x=0.03 in Prx_2.La.jPb3. 12j NFL behavior appearing in the specific heat for x>0.95 
was attributed to the quadrupole Kondo effect. The absence of AFQ order between the T 3 
doublets in PrlnAgo and PrMg3 was discussed by taking the quadrupole Kondo effect into 



consideration. 
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We have recently reported that a cubic compound PrIr 2 Zn 2 o undergoes an AFQ order 



ing at 7SHM1 K and a superconducting transition at 2-0.05 K.|l_ PjThis compound 
crystallizes in the cubic CeCr 2 Al 2 o-type structure with the space group Fd3m and Z=8.[19j 
where the Pr atom is encapsulated in a highly symmetric Frank-Kasper cage formed by 
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16 zinc atoms. The superconductivity below T c was indicated by the sizable diamagnetic 
signal in the AC magnetic susceptibility. Our analysis of the magnetic anisotropy in the 
paramagnetic state and a Schottky peak in the specific heat revealed that the CEF ground 
state is the non-Kramers T 3 doublet and the two-fold degeneracy is released by the AFQ 
ordering at Tq=0.11 K. Thus, this compound is the first example where the superconducting 
transition occurs in the AFQ ordered state. The entropy at Tq is reduced to 0.2i?ln2 from 
i?ln2 that is expected for two-fold degeneracy of the r 3 state. This strongly suggests that 
the quadrupole fluctuations play a role in forming the superconducting Cooper pairs. 

Matsushita et al. have reported that the specific heat has no peak at Tq but a broad 
peak at around 0.4 K for their single- crystalline sample of PrIr 2 Zn 2 o 20]. It was interpreted 
as the manifestation of a quadrupole-glass-like ordering. The absence of the sharp peak due 
to quadrupolar transition, however, might result from the lower quality of their sample than 
that of ours. In fact, the residual resistivity of their sample is 40 times larger than ours. 
From the very small effective mass at the Fermi level, they argued very weak hybridization 
effect of the 4/ electrons with the conduction electrons. 20] 

It is worthy noting that for our samples of PrIr2Zn2o the specific heat divided by tem- 
perature, C/T, shows the — InT dependence and the electrical resistivity p(T) shows VT 
dependence in the range between 0.2 K and 0.8 K. 23jj These temperature dependences are in 
accord with the theoretical calculation based on a two-channel Anderson lattice model with 
a low characteristic temperature. |22J This good accordance suggests that quadrupolar exci- 
tations are weakly coupled to the conduction electrons. We note that the above temperature 
dependences in C/T and p(T) have been reported in an isostructural compound PrV^A^o 
with the T 3 doublet CEF ground state above Tq. 23[ Furthermore, another isostructural 
compound PrTi2Al2o has been found to show a superconducting transition at 0.2 K in the 
ferroquadrupole ordered state. Thereby, clean limit superconductivity and enhanced effec- 
tive mass of m*/m~ 16 have been suggested. 24J Therefore, in order to address the issue 
whether the superconducting pair is mediated by the quadrupole fluctuations or not, further 
experimental investigations of magnetic and transport properties of the family of PrT 2 X2o 
(X=A1 and Zn) are certainly needed. 

In the present paper, we focus on PrRh 2 Zn 2 o which is isoelectronic to PrIr 2 Zn 2 o- Be- 
cause similar CEF level scheme is expected for the two systems, the CEF ground state in 
PrRh 2 Zn 2 o would be the T 3 doublet. It is intriguing to know whether the AFQ ordering 
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and superconducting transition occur or not in PrRli2Zn2o. Keeping this in mind, we have 
grown single- crystalline samples and studied the magnetic and transport properties. A part 
of the present work has been published in conference proceedings. 2l|, |25] 



II. EXPERIMENTAL 

Single-crystalline samples of PrRh2Zn2o were grown by the zinc self- flux method with 
high purity elements of Pr (4N), Rh (3N) and Zn (6N). The CeC^A^o-type structure was 
confirmed by the powder x-ray diffraction analysis. By X-ray diffraction analysis for a 
single-crystalline sample at 293 K with Mo Ka radiation using a Bruker SMART APEX 
CCD area-detector diffractometer, the lattice parameter at room temperature was refined 
to be 14.2702(3) A, whose value is a little smaller than 14.287(1) A in the previous report. [2] 



It is comparable with that for the isoelectronic PrIr2Zn2o, 14.2729(2) A. 16] The chemical 
compositions of the crystals were determined by averaging over 10 different regions for each 
crystal with a JEOL JXA-8200 analyzer. The compositions of Pr : Rh : Zn = 1 : 1.98(1) : 
20.02(5), where the numbers in the parentheses are the standard deviations. The electron- 
probe microanalysis revealed the impurity phases of binary alloys RhZn6.4 and Pr^Zn whose 
volume fractions are less than 5% in view of back-scattering electron image. The single- 
crystalline sample was oriented by the back reflection Laue method using an imaging plate 
camera, IPXC/B (TRY-SE). 

The electrical resistivity p was measured by a standard four-probe AC method in a 
laboratory-built system with a Gifford-McMahon-type refrigerator between 3 and 600 K. 
The measurements down to 0.045 K were done in a laboratory-built system with a commer- 
cial Cambridge Magnetic Refrigeration mFridge mF-ADR50. Magnetization was measured 
using a commercial SQUID magnetometer (Quantum Design MPMS) between 1.9 and 300 
K in magnetic fields up to 5 T. The AC magnetic susceptibility was measured down to 
0.045 K using a laboratory-built system installed in the mFridge mF-ADR50. The mea- 
surements were carried out for a frequency of 1 kHz, with an applied AC field of 0.02 mT. 
Single- crystalline samples of PrRh 2 Zn 2 o and the reference superconductor LaIr 2 Zn 2 o of the 
same dimensions 0.5x1x4 mm 3 were loaded into pickup coils. The specific heat C was 
measured by a relaxation method using a commercial calorimeter (Quantum Design PPMS) 
at temperatures between 0.4 and 300 K. The measurements down to 0.047 K were done by a 
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quasi-adiabatic method with a 3 He- 4 He dilution refrigerator. Thereby, horizontal magnetic 
field up to 7 T was applied by a split-pair superconducting magnet. In zero magnetic field 
measurements, the dilution refrigerator was inserted into a cryo-dewar without a supercon- 
ducting magnet to avoid residual magnetic fields. 



III. RESULTS AND DISCUSSION 



Figure [T] shows the temperature dependence of the electrical resistivity p(T) of PrRhjZ^o 
between 30 K and 600 K. The hysteretic behavior appears in a wide range between 140 and 
470 K, which indicates a first-order phase transition. However, the specific heat data shown 
in the inset of Fig. [1] exhibits no apparent peak at around T s , because the first-order phase 
transition is hard to be observed by the relaxation method. Hysteretic behaviors in p(T) 
were also observed in LaRu2Zn 20 and PrRu2Zn 2 o at T s =150 and 138 K, respectively. 16J 
Since superlattice reflections were observed in PrRu2Zn2o below T s by the electron diffrac- 
tion method, the hysteretic behavior was attributed to a structural transition. 16J On the 
other hand, no evidence for structural transition was found in p(T) from 530 K to 30 
K for PrlroZ^n- Therefore, the quadrupolar degrees of freedom remain active in the 
doublet. 17[ It undergoes antiferroquadrupole ordering at Tq=0.11 K and superconducting 



transition at T c =0.05 K. 

Figure [2] shows the low-temperature part of p(T) of PrRh2Zn2o in various constant mag- 
netic fields .8=0, 1.5, and 5.0 mT applied along the [111] direction. In B=0, p drops to zero 
at T c =0.06 K, indicating a superconducting transition. The bulk nature of the superconduc- 
tivity was confirmed by detecting a large diamagnetic signal of the Meissner effect below T c 
in the AC magnetic susceptibility as shown in the inset of Fig. [2j The diamagnetic signal 
below T c for PrRh 2 Zn 2 o is about half that of the reference superconductor LaIr 2 Zn 2 o with 
T c =0.6 K, indicating that the superconducting volume fraction in the sample is about 50 
% at the lowest temperature of 0.045 K. Applying magnetic fields of 1.5 mT and 5.0 mT, 
the drop of p is suppressed and the normal state remains to the lowest temperature 0.05 K. 
This means that T c is suppressed below 0.05 K even in the small magnetic field of B=1.5 
mT. The vanishment of T c in B—1.5 mT seems to be consistent with the critical magnetic 
field of 0.5 mT predicted from the BCS theory. The magnetic field dependence of T c at 
low fields below 1 mT should be measured in detail to evaluate the critical field. Below 
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0.4 K, the p(T) increases with increasing magnetic fields, probably because fluctuations of 
field-induced magnetic moments effectively scatter the conduction electrons. 

The data of x{T) for PrRh 2 Zn 20 is plotted in the lower inset of Fig. El On cooling below 
10 K, x{T) gradually approaches a constant value, indicating a Van-Vleck paramagnetic 
ground state. Thus, the CEF ground state is nonmagnetic, either F 1 singlet or T 3 doublet. 
The main panel of Fig. [3] shows x _1 (T) measured at 5=0.1 T. Above 30 K, follows the 
Curie- Weiss law with the effective magnetic moment of 3.55 /ie/f-u., which is in agreement 
with the value of the trivalent Pr free ion. The small value of —4.8 K for the paramagnetic 
Curie temperature 9 P indicates that the intersite magnetic interaction between Pr ions is 
antiferromagnetic but rather weak. 

To judge whether the CEF ground state is the nonmagnetic either Ti singlet or T 3 doublet, 
we measured the isothermal magnetization M(B) at T—1.8 K by applying magnetic fields 
up to 5 T along the [100], [110], and [111] directions. As shown in the upper inset of Fig. 
El the M(5)'s data are almost the same up to 2 T, above which three curves gradually 
diverge; M(B || [100]) > M(B || [110]) > M(B || [111]). This anisotropic behavior can be 
well reproduced by the calculation of M(B) using the r 3 — CEF model, as shown with 
the solid lines. The T 3 — T 4 level scheme separated by 32 K is depicted in the inset. Figure 
H] shows the temperature dependence of the magnetic part of the specific heat divided by 
temperature, C m /T. For the lattice part of C in PrRh 2 Zn 20 , we used the data of LaIr 2 Zn 20 
because the more relevant compound LaRh 2 Zn 2 o does not form. As shown in the inset 
of Fig. HI a broad peak manifests itself at around 10 K. To reproduce this Schottky-like 
anomaly, we used the T 3 — model. The solid line represents the calculation, where 
the higher excited levels were not taken into account. The good fit corroborates the T 3 
doublet ground state. The T 3 ground state has been also supported by the inelastic neutron 
scattering experiments, where the excitations from the T 3 ground state to the excited states 



were observed. 



261 ] Furthermore, no additional excitation peak appears below T s , indicating 



that the cubic point group of the Pr site is conserved even below T s . 

The gradual release of the entropy of the T 3 doublet results in the gradual increase in 
C m /T for T<3 K. Between 0.6 and 0.1 K, C m /T obeys — InT dependence as is shown with 
the solid line in Fig. [4j It is probably the manifestation of the two-channel Kondo effect 



due to the quadrupolar degrees o 



PrT 2 Al 20 (T=V and Nb).|21 



,l23j 



freedom as were found in the isostructural PrIr 2 Zn 20 and 
The slope of the —InT dependence of C m /T in Fig. @] 
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gives the characteristic temperature Tk as 0.7 K.j4] The entropy release at Tk from i?ln2 is 
1.35 J / K 2 mol, which value is moderately consistent with the theoretical value 1.45 J / 
K 2 mol calculated by taking the two-channel Kondo effect into consideration.^] On further 
cooling, a peak appears at 0.06 K, which coincides with the superconducting transition at 
T c =0.06 K. To judge whether this peak in C m /T is due to the superconducting transition 
or not, we measured the specific heat in magnetic fields. The inset of Fig. [5] (b) shows 
the data in B=0, 0.3, and 0.5 T applied along the [110] direction. The peak at 0.06 K in 
B=0 slightly shifts to higher temperatures at 5=0.3 and 0.5 T, indicating the phase tran- 
sition remains even in B=0.5 T. Note that the resistive transition into the superconducting 
state is quenched by the very weak field of 1.5 mT as shown in Fig. [2J Assuming a BCS 
superconductor with electronic specific heat 7= 10 mJ/K 2 mol, the jump in the specific 
heat, AC, is evaluated as 8 mJ/K mol at T c =0.06 K from the relation of AC/7T C =1.43. 
This value of AC is smaller than the resolution of our specific heat measurements with the 
quasi-adiabatic method, say 50 mJ/K mol at 0.06 K. Therefore, we conclude that the peak 
of the specific heat at 0.06 K does not originate from the superconducting transition but a 
phase transition due to the multipolar degrees of freedom in the T3 doublet ground state. 
Thus, we mark this transition temperature as Tq hereafter. The magnetic entropy S was 
evaluated from the temperature dependence of the C m /T as shown with the solid curve in 
Fig. HI On cooling, the S(T) monotonically decreases from the value of i?ln2 at around 2 
K. The S at Tq=0.06 K is only 0.1i?ln2 which is much smaller than i?ln2 that is expected 
for the two- fold degeneracy of the T 3 doublet. 

Figure |5] shows the temperature dependence of the specific heat in magnetic fields applied 
along (a) [100], (b) [110], and (c) [111]. For B || [100], Tq decreases with increasing the 
magnetic field, and disappears above 3 T. On the other hand, for B || [110] and [111], Tq's 
increase with increasing the magnetic fields up to 6 T. Above 5 T, Tq starts decreasing 
in B || [110], whereas Tq still increases in B || [111]. The magnetic field dependences of 
Tq are summarized in the B—T phase diagrams in Fig. |6] (a), where AFQ, SC, and PM 
mean the antiferroquadrupolar ordering, superconducting, and paramagnetic states, respec- 
tively. In B=0, the superconducting transition and AFQ transition occur simultaneously at 
T c =Tq=0.06 K. Note that Tq depends on the magnetic field directions; Tq (B || [100]) < 
Tq (B II [110]) < Tq (B II [111]). This relation is similar to those observed in PrPb 3 and 
PrIr2Zn2o with the non-Kramers T3 doublet ground states. 
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The anisotropic behavior of Tq(B) in PrRh2Zn2o can be explained by the calculation of 
the CEF level scheme in the paramagnetic state. The splitting of the r 3 doublet calculated 
for B || [100], Aioo, is much larger than A 110 in B 1 1 [ 1 1 0] and A m in B ||[111]. This 
calculation reproduces the isothermal magnetization at 1.8 K shown in the upper inset of 
Fig. El where M(B || [100]) exceeds M(B || [110]) and M(B || [111]) above 2 T. If the 
energy scale of the quadrupole interaction becomes smaller than the field induced splitting 
of the r 3 doublet, the quadrupole freedom would not order any more. Therefore, the order 
of A ioo > A 110 > A m in the magnetic fields corresponds to the order of Tq; Tq (B || [100]) 
<T Q (B I) [HO]) < Tq (B I) [111]). 

Furthermore, for the quantitative analysis of the above results, we performed mean-field 
calculation with a two-sublattice model based on the following isotropic interactions: 

"^A(B) — ^CEF - QjUbJH - (A'i(J)b(A) + K 2 {J) A(B)) J 

- K r3 [(0° 2 ) B{A) 0° 2 + (0 2 2 ) B( A)0 2 2 ], 

(1) 

where gj is the Lande ^-factor (= 4/5 for Pr 3+ ) and //b is the Bohr magneton. The first 
term in eq. ([T]) is the CEF Hamiltonian for the subspace of the J=4 multiplet, which can 
be written as 

no _ c;n 4 no _ oi n 4 

^ = Wb%^ + {l-W\)^o], (2) 



60 

where the notation by Lea et al. is used. 



2n 



Here, we adopt the two parameters W=— 1.1 



K and £=0.46 which were determined to reproduce both the magnetic susceptibi 



ity as 



shown with the solid curves in Fig. 3 and the inelastic neutron scattering spectra[26j. The 
eq. © leads to the CEF level scheme as T 3 (0) - T 4 (31 K) - T 5 (65 K) - T t (73 K). 
In eq. ([1]), K\ and K 2 are the inter- and intra-sublattice magnetic interaction coefficients 
for excited magnetic multiplets, respectively, and Kr 3 is the inter-sublattice interaction 
coefficient of r 3 -type quadrupolar moments. It should be noted that the quadrupoles 2 
and 0\ have the same interaction coefficient in this model. As a result, the 2 and the 0\ 
phases are degenerate for B — 0, i. e. Tq(0 2 ) = Tq(0 2 ). From the experimental result 
Tq=0.06 K for B=0, we obtain K-p 3 =— 0.0037 K. The difference between the measured 
magnetic susceptibility x±f a ^ T>Tq and Xcef calculated using the CEF parameters shown 
above leads to the mean-feild parameter A=— 1.5 mol/emu, which represents the exchange 
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interaction among the Pr magnetic moments, and the relation K\+K2=— 0.4 K. Although 
these antiferromagnetic interactions are much stronger than the AFQ interaction, the system 
does not order antiferromagnetically because the non-magnetic T 3 doublet is well separated 
from the excited magnetic multiplets. 

The calculated B—T phase diagrams for B || [100], B || [110], and B || [111] are shown 
with the thick curves in Fig. [6] (b). For simplicity, we set -fT 2 =0 in the relation K1+K2—— 
0.4 K. For B || [100], [110], and [111], Tq increases with increasing magnetic fields up to 
B=l, 3, and 5 T, and at T=0 the AFQ ordered phase closes at the critical fields of B=2.3, 
4.0, and 6.3 T , respectively. This calculation of Tq(B) qualitatively reproduces both the 
anisotropic response of Tq and the initial increment of Tq with the magnetic fields. Although 
the calculated critical fields of the AFQ ordered phase are much lower than the experimental 
data. The quantitative disagreement is owing to simple two-sublattice model where only the 
isotropic magnetic and quadrupole interactions were taken into consideration. Quadrupole 
structures could be more complicated than the two-sublattice model and there exist various 
kinds of interactions between the higher-order multipoles induced by magnetic field. In 
order to deduce the effect of the antiferromagnetic interaction on Tq(B), we calculated Tq 
by setting Ki=K 2 =0. As is shown by the dashed curves in Fig. [6] (b), Tq's do not increase 
with the magnetic field and the AFQ ordered phase collapses at rather low fields of 1.8, 2.7, 
and 4.1 T for B || [100], [110], and [111], respectively. This discrepancy strongly suggests 
that the AFQ ordered phase is stabilized by the AFM interaction between the field-induced 
magnetic dipoles. 

The important fact of the B—T phase diagram of Fig. |6] (a) is that the superconducting 
region is inside the AFQ ordered phase as found for the isoelectronic PrIr 2 Zn 2 o. jl3] As shown 
in Fig. HI the S(T) at Tq is reduced to 0.1i?ln2 from i?ln2 that is expected for the two- 
fold degeneracy of the T 3 state. The reduced entropy at Tq in PrRh 2 Zn 20 from 0.2i?ln2 in 
PrIr 2 Zn 20 is attributed to the quadrupole fluctuations. Although T c =0.05 K is lower than 
Tq=0.11 K in PrIr 2 Zn 2 o, [itJ the two transitions simultaneously occur at T c =Tq=0.06 K in 
PrRh 2 Zn 2 o- It is interesting to examine if the quadrupole fluctuations would stabilize the 
superconducting state. To reveal how the superconductivity couples with the quadrupolar 
degrees of freedom, a systematic study of the family of PrT 2 X 20 (X=A1 and Zn) at low 
temperatures and under high pressures are needed. 

It remains as an important issue why the CEF ground state in PrRh 2 Zn 2 o is the r 3 dou- 
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blet state in spite of the possible structural transition between 140 and 470 K. Because the 
T3 doublet is realized only in the cubic point groups, the low temperature phase must belong 
to another cubic point group. We recall that the metal-insulator transition in PrRu4Pi2 at 
Tui— 63 K is accompanied with a structural transition from a body-centered cubic struc- 
ture in the high-temperature metallic phase to a simple cubic one in the low-temperature 



insulating phase [29j. In order to determine the structure of the low-temperature phase 
of PrRli2Zn2o, x-ray diffraction and electron diffraction measurements on single-crystalline 
samples are in progress. 



IV. CONCLUSION 



We performed electrical resistivity p, magnetic susceptibility Xi an d specific heat C mea- 
surements on PrRh 2 Zn 2 o, which is isoelectronic to PrIr 2 Zn 20 showing antiferroquadrupolar 
and superconducting transitions at Tq=0.11 K and T c =0.05 K, respectively. In PrRh 2 Zn 2 o, 
we found a superconducting transition at T c =0.06 K, below which the bulk nature of the 
superconductivity was confirmed by the large diamagnetic signal due to the Meissner effect. 
The analysis of both the Schottky anomaly in specific heat and anisotropic magnetization 
curves for B>2 T at T—1.8 K indicates that the CEF ground state of the trivalent Pr ion 
is the non-Kramers T 3 doublet with the quadrupolar degrees of freedom. The sharp peak in 
the specific heat at 0.06 K is attributed to the antiferroquadrupolar ordering. This temper- 
ature Tq shows anisotropic response to magnetic field B; Tq decreases in B || [100], whereas 
Tq increases in B || [110] and B || [111] with increasing B up to 6 T. This anisotropic 
behavior of Tq(B) can be well explained by the mean- field calculation with a two-sublattice 
model, assuming the isotropic inter-sublattice interaction between the T3-type quadrupolar 
moments and the magnetic interaction between the field-induced magnetic dipole moments. 
These experimental and calculated results corroborate the AFQ ordering below Tq. The 
superconducting transition and AFQ order occur simultaneously at T c =Tq=0.06 K in B=0. 
By applying magnetic fields, however, the superconducting state is easily quenched at several 
mT but the antiferroquadrupolar phase is stabilized. To answer the issue whether the su- 
perconducting pair is mediated by the quadrupole fluctuations or not, further experimental 
and theoretical works are needed. 
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FIG. 1: (Color online) Temperature dependence of the electrical resistivity of PrRli2Zn2o in the 
temperature range between 30 K and 600 K. Tg own , Taf and denote the temperatures of 
anomalies in the cooling and heating processes, respectively. The inset shows that the specific heat 
does not exhibit any anomaly at T$. 
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FIG. 2: (Color online) Temperature dependence of the electrical resistivity p in B=0, and in 
magnetic field of 1.5 and 5 mT applied along the [111] direction. The inset shows the AC magnetic 
susceptibility of PrRli2Zn2o and the reference compound LaIr2Zn2o with T c =0.6 K. 
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FIG. 3: (Color online) Temperature dependence of the inverse magnetic susceptibility \~ l °f 
PrRli2Zn2o- The solid curve is a fit of the data with the CEF parameters of W=— 1.1 K and 
a;=0.46, and the mean-field parameter A=— 1.5 mol/emu. The lower inset shows the magnetic 
susceptibility x(^)- The upper inset shows the isothermal magnetization at 1.8 K in the magnetic 
fields applied along the [100], [110], and [111] axes. 
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FIG. 4: (Color online) Temperature dependence of the magnetic specific heat divided by temper- 
ature, C m /T (left-hand scale) and the entropy (right-hand scale). The inset shows the C m /T on a 
linear temperature scale. 
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FIG. 5: (Color online) Temperature dependence of the specific heat in various constant magnetic 
fields for (a) B || [100], (b) B \\ [110], and (c) B || [111]. 
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FIG. 6: (Color online) (a) B — T phase diagrams for B || [100], [110], and [111] determined by the 
specific heat measurements. Tq and T c are the AFQ and superconducting transition temperatures, 
respectively. AFQ, SC, and PM indicate the antiferroquadrupolar ordering, superconducting, and 
paramagnetic states, respectively, (b) B — T phase diagrams obtained by the two-sublattice mean- 
field calculation. The CEF parameters are W=— 1.1 K and x=0.46, and the quadrupole interaction 
is Kr 3 =— 0.0037 K. The solid and dashed curves indicate the phase boundaries calculated with the 
intersite magnetic interaction K\=— 0.4 K and without the magnetic interaction Ki=0, respectively. 
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